Wave-packet analysis of interference patterns in output coupled atoms 
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We study the output coupling of atoms from a magnetic trap into a linear potential slope of 
gravity using a weak radio-frequency field. We present a one-dimensional wave-packet model based 
on a continuous loading of a continuous spectrum of generalised eigenstates to describe the scenario. 
Analyzing the model, we show how the interference of the classical coupling fields maps to the 
interference of the resulting atomic streams. 
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I. INTRODUCTION 



ence patterns due to multiple simultaneous couplings. 
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Ever since the first realization of an atomic Bose- 
Einstein condensation [jJ-Q, there have been applica- 
tions, where the coherent cloud of trapped atoms has 
been used as a source for output coupling Q . The coher- 
ence properties of the source can be mapped to a coherent 
output and, moreover, by applying continuous cou- 
pling a coherent stream of spatially wide-spread atoms 
can be created 0, [l(| ■ In a close analogy with the opti- 
cal laser, an atom laser is thus formed. 

An often used method to realise such a system, is to 
induce spin flips to the trapped cloud of atoms by in- 
troducing a weak magnetic field perturbation, i.e., an 
oscillating radio frequency (rf) field perpendicular to the 
static trapping field [H, |Tl[ H^ . The rf- field creates a 
coupling between the Zeeman sublevels Mp and the in- 
ternal spin state can be thus nipped to an untrapped 
or, with strong rf-field intensities, even to anti-trapped 
states [HI, [T3 |. Especially, in the linear Zeeman shift 
regime, the sublevel Mp = does not couple to the static 
trapping magnetic field at all, but is only affected by the 
linear potential slope of the gravity. Consequently, such 
atoms fall freely and exit the trapping area. Other im- 
plementations of the output coupling are, e.g., to apply a 
Raman transition [l5[ , which could also provide the free- 
falling atomic flux with an initial momentum kick, or to 
construct a tunneling connection [l||. Interestingly, the 
output coupling situation is reminiscent of the molecular 
dissociation triggered by ultrashort pulses [l7l - [l9j . 

Over the years, there have been theoretical papers con- 
sidering the atom lasers with one-dimensional (20l - [22l | 
and three-dimensional model s [l2l. [231 |24| . using weak 
(HIMIIfH and strong [HlSllillJcouplingstrengths, 
applying multiple simultaneous couplings (2lll23j , having 
the source at finite temperature [25|, [26| , and from the 
point of view of stability [2714291 ] and pulse shape [3(| • In 
this paper we analyze a simple one-dimensional model in 
order to clarify one specific problem concering interfer- 
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The phase coherence of the spatially elongated atomic 
beams is most strikingly demonstrated by strong interfer- 
ence patterns while superimposing two beams with dif- 
ferent energies [3ll - l33j . Again in a close analogy with the 
optical lasers, the interference pattern depends on three 
quantities: (i) the relative amplitudes, (ii) the relative 
phase difference, and (iii) the energy separation. 



The spatially wide-spread wave- function interpretation 
of the interfering atomic beams is a strongly non-classical 
result. However, an alternative explanation in terms of 
interfering (classical) magnetic rf-fields, which are driving 
the coupling, has been proposed 33] . In this line of rea- 
soning the coupling magnetic field is understood in terms 
of a carrier frequency and a beating envelope, and the 
correspondence between the pulsing rf-amplitude and the 
resulting output stream was demonstrated. There seems 
to be a discrepancy between the two ways of looking at 
the problem. On the one hand, the system is described by 
a pulsing flux generated by a pulsing semi-classical cou- 
pling, while on the other hand, the system is described by 
interference of superimposed spatially elongated asymp- 
totic atomic wavefunctions |2ll |3l| . The purpose of this 
paper is to demonstrate the connection between these 
two extreme interpretations. 



In Sec. [II] we derive a wave-packet solution to a sim- 
plified one-dimensional problem in terms of a continuous 
loading of a continuous spectrum of generalised energy 
eigenstates. In Sec. IIIII we show how the visibility of 
the atomic interference pattern maps from the interfer- 
ence of the magnetic fields. We then apply the model in 
Sec. IIVI using realistic experimental parameters and com- 
pare the results with numerical simulations including the 
complete Zeeman-sublevel structure as well as the atomic 
contact interactions. Finally, we finish with conclusions 
and discussion in Sec. [V] 
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II. WAVE PACKET MODEL 

A. Physical system 

An atom couples to the magnetic field via its magnetic 
moment, resulting in an interaction energy defined as 



[/(B) = -/x-B. 



(1) 



Above, the magnetic moment operator is fi = —fio(gsS + 
giX> + 9iI)/K where fio = \e\h/2m e is the Bohr mag- 
neton and gi are the Land g-factors for electronic spin 
(S), orbital (L), and nuclear spin (I) angular momentum. 
When the energy splitting corresponding to this term is 
small compared to fine and hyperfine splittings, the total 
angular momentum F = I + J, with J = L + S, is a good 
quantum number and fi ~ —fiog F F/h, where the Land 
factor is g F ~ gj[F(F+l)+ J(J+l)-7(/+l)]/2F(F+l), 
with gj ~ 1+[J(J+1)+S(S+1)-L(L+1)]/2J(J+1). In 
the limit of weak magnetic field, which is the case in the 
present work, the Zeeman splitting between the sublevels 
M F is linear fjl . 



with a finite component in the direction perpendicular 
to the trapping field. The pulse envelope B (t) has an 
arbitrary shape and the pulse is turned on after the initial 
time t = 0. As will be clear from the following, the model 
can be generalised directly to any linear combination of 
such single-mode rf-fields. Consequently, it is sufficient 
now to consider a single rf-field. 

The rf-field results in an interaction Hamiltonian 



H T (t) = -n ■ B rf (t). 



(4) 



We will write the polarization vector as e r f = 
52i=+ t - <z (&i ■ e r f)e i; where e± = {e x ± ie y )/s/2. The z- 
component causes only a small perturbation to the trap- 
ping potential, and is assumed to be zero hereafter. The 
circular components, corresponding to the raising and 
lowering angular momentum operators F± = F x ± iF y , 
induce transitions between the sublevels, as F± \F, M F ) = 
H^F(F + 1) - M F (M F ± l)\F, M F ± I). Finally, we re- 
mark that the total angular momentum operator F 2 com- 
mutes with the total Hamiltonian H(t) = Hq + Hi(t) as 
well as its components Hq and Hi(t), so the dynamics is 
confined into a single hyperfine state F. 



1. Trapping potential and gravity 



B. Representation of the state 



A magnetic trap for the atoms in the low-field-seeking 
states is formed by simply creating a magnetic field inten- 
sity minimum. The local direction of the field describes 
the quantization axis e z , and close to the minimum, the 
magnetic field is assumed to be approximately harmonic, 
such that B trap = Bt ra p (r)e z = B^. ap (\lx 2 + \ 2 y y 2 + 
\ 2 z 2 )e z . In the same direction, a strong static bias field 
Bbias = -BbiasGz is applied in order to remove the degen- 
eracy at origin and, hence, to suppress the Majorana spin 
flips and the resulting atom losses 5] . 

The trapping potential operator is C/trap(r) = 
sga{g F )[\m{uj 2 x x 2 + uj 2 y 2 + uj 2 z z 2 ) + hw hias ]F z /h, where 
•J 2 = 2/j, \g F \B^ ap \ 2 /m and u hias = fi Q \g F \B hias . The 
atoms feel, irrespective of their internal state, also the 
linear potential of the gravity, C/ grav ity(r) = —mgx; the 
harmonic trapping potentials are relocated accordingly 
in position and energy. The static Hamiltonian reads 



H Q = T + t/trapO) + ^gravity (r), 



(2) 



where T = — h 2 V 2 /2m is the kinetic energy term. As 
is now obvious, an integer- valued hyperfine state F sup- 
ports a special sublevel M F — 0, which feels only the 
linear gravitational potential. 



2. Coupling rf-field 

The coupling between the Zeeman sublevels is induced 
by applying a weak rf-field 
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The coupling is assumed to be weak, so only transitions 
to sublevels M^gnai = M F .i n iti&\ ± 1 are relevant. Since 
the magnetic (trapping) potentials for the different sub- 
levels are (F, Mi?|C/tra P |F, M F ) oc sgn(g F )M F , we will as- 
sume that the trapped atoms are initially on the internal 
state |T) = \F,M F = sgn(g F )) and that the free-falling 
untrapped state is \U) = \F,M F = 0). The transition 
between the internal states |T) — » \U) is provided by the 
operator F + /_ in systems with negative/positive g F . 

In the following, we will neglect any population trans- 
fer to the antitrapped high-field-seeking Zeeman sub- 
levels [M F = —nsgn(g F ), with n — 1, . . . , F], for which 
the magnetic field minimum forms a repulsive poten- 
tial, as well as to the more energetic trapped sublevels 
[M F = nsgn(g F ), with n = 2, . . . ,F]; the relevant po- 
tentials and couplings are illustrated in Fig. [TJ This is 
a justified neglection since we are interested in the weak 
coupling regime. On the other hand, if a maximal flux of 
atoms would be desirable, one would have to use strong 
rf-fields, and in that case these neglected sublevels would 
have a nontrivial contribution LJj. With ever stronger 
coupling strengths th e sy stem should be described by 
dressed potentials |35l - [38| |. 

The convenient choice of basis functions depends on 
the internal state. For the trapped state |T), the ba- 
sis is provided by the harmonic oscillator eigenstates 
{ | (/>„)} ^ =0 . For the untrapped state \U), however, the 
external potential is linear, and the basis is formed by 
an uncountable set of generalised eigenfunctions, or dis- 
tributions, \\iI> f )}epr, which satisfy the Airy differential 
equation [391 ] . Their explicit form is 



B r f(i) = -B Q (t)e l{ e 



-i{u Ti t+9) 



+ C.C. 



(3) 



ijj E (x) =AfAi[(x + E/mg)/l], 



(5) 
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FIG. 1: (Color online) Schematic setup of potentials and cou- 
plings. The trapped Zeeman sublevel \T) is coupled to the un- 
trapped level \U) by a weak rf-field. Starting from a trapped 
state with energy Eq, the rf- frequency determines the reso- 
nance energy -E rcs , around which a continuous spectrum of 
generalised energy eigenstates is populated during the output 
coupling, (a) With multiple simultaneous frequency compo- 
nents Wrfj, a corresponding set of resonant energy levels is 
formed. The atomic streams interfere as they fall in the grav- 
ity field, (b) Equivalently, the coupling can be interpreted as 
being driven by the sum of the single rf-fields, which corre- 
spond to a carrier frequency cj r f,+ = (w r f,i+w r f,2)/2 and a cou- 
pling strength pulsing at frequency cj r f,- = (wrf.i — <^rf,2)/2. 

where the normalisation factor J\f = X/ly/mg and the 
characteristic length scale I — (H 2 /2gm 2 ) 1 ' 3 . These gen- 
eralised functions are not normalisable according to the 
L 2 norm, (iPe\?Pe'} = S(E — E'), and thus can not in- 
dividually represent any physical state. However, they 
form a complete orthonormal spatial basis, in a sense 
that / &E ip* E {x)ip e{x') — S(x — x'). Therefore, any spa- 
tial state \<p) can be described in terms of these distri- 
butions as \<p) = J dE f (E)\ip e) , where the spectrum is 
f{E) = (i/je\(p). Consequently, the normalisation of the 
state is done in accordance with the properties of the 
spectrum, such that \\tp\\ 2 = J dE\f(E)\ 2 . It is imme- 
diately evident, that any state \ip£>,U) described by a 
discrete spectrum f(E) = J2 { CiS(E — £7,-) is, first of all, 
unphysical, and corresponds to a (quasi)periodic solution 
within the evolution generated by Ho- This does not fit 
to the intuition about a free-fall event. Combining the 
previous statements, any state within our system can be 
expressed as 

!*(*)> =£M*)l&» r >+ / dEc B (t)\i> B ,U). (6) 

In the interaction picture with respect to Hq, given 
by Eq. ©, = e iHD */ ft |*(i)) and the correspond- 

ing coefficients are b n (t) — e lEnt l h b n it) and Cs(t) = 
e lEt l h CE{t\ In the following, we will also use the no- 
tation 1*0) = where /3 = T,U, for the trapped 
and untrapped components of the total state. 



C. Wave-packet solution 

The coupling between the different sublevels comes 
from the interaction Hamiltonian Q. The weak cou- 
pling causes only a small perturbation to the bare sys- 
tem, defined by the static Hamiltonian ^ , and therefore 
its effect can be described by transition matrix elements. 
Let us assume that initially the system is at equilibrium 
in a trapped ground state |^(0)) = \<j>,T), for which 
Ho\<f>) = Eo\(f>). In terms of the representation ([5]), the 
coefficients are &o(0) = 1, b n (0) — c\e(0) = for all n > 
and £el. 

In the interaction picture the equation of motion for 
the coefficients cg(i) is given by 

P E (t) ^EM^{t))^- l -^ E ,U\Hi{t)\^{t)) 

= - '-bo^e-^-^^El^iUlHj^lT). 
n 

(7) 

The trapped state remains essentially intact during the 
weak coupling pulse, so we can assume bo (t) = 1. The 
formal solution in the Schrodinger picture is 

c E {t)=-U- iEt ' h (M<t>) f dse-^- E >/ h 
ft Jo 

x (U\Hj( S )\T). (8) 

Therefore, according to definition ([5]), the untrapped 
component is given by 

\^u{t))=-\f^se-^ s ' h (U\H I { S )\T) 

x J dEe-^-^&Em^E). (9) 

Defining the outcoupling rate function Q and the respec- 
tive instantaneous outcoupled state |$), corresponding 
to a delta-peak outcoupling rate function, as 

n(t) = ~le- lEnt / h (U\H I (t)\T), (10) 
n 

|*(*)) = J dEe- iEt / h &E\<P)\i>E) = (Ule-^/^U), 

(11) 

the full time-dependent solution for the outcoupled 
atomic beam can be written in a compact form as a con- 
volution 

|*et(*)> = / d8n(8)\$(]t-8))= [n*(0|$))](t), (12) 
Jo 

where the Heaviside theta function, for which 0(t) equals 
to for t < and 1 for t > 0, takes care of a proper 
temporal causality. The instantaneous outcoupled state 
|3>(£)) matches the static-Hamiltonian-induced evolution 
[cf. Eq. @] of the spatial component of the initial 
trapped state | <$>) , only its internal state is the untrapped 
one. Finally, we remind that f2(£) vanishes for t < 
according to our previous definition. 
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D. Continuous spectrum of states 



Interference of classical fields 



Let us consider the matrix element of the interaction 
Hamiltonian between the trapped and untrapped states 
(U\Hj\T). Since the trapped and the untrapped states 
are separated by a single quantum of angular momen- 
tum, (T\F Z \T) = sgn(g F )h = ±H and (U\F Z \U) = 0, the 
transition between the states \T) — > \U) is induced by 
the operator F a , where a = +/— for systems with nega- 
tive/positive gF- Therefore, Eq. ([7]) is 

j t c E (t) = - ^-J^M9FB (t)(rp E \cf>)(U\F a \T) 
x [(e Q • e T{ ) e ^KEo+^-E)t+e] 

+ (e a • g* f ) e -«[(Eo-" rf -B)t-e]j ( (13) 

where the factor \J~2 comes from the identity e± • F = 
F±/V2. 

With a constant rf- field, Bo(t) = Bo<d(t), the time 
integration gives terms 



c E (t) oc t{ipi 



. {E Q ±LO Ti ~E 
>sinc 1 



(14) 



Therefore, the spectrum concentrates in the vicinity of 
resonant energy levels E = Eq ± w r f as time passes. 
According to the physical setup, on the other hand, 
the overlap integral (ipE\<fi) is concentrated around E ~ 
— mg{4>\x\(j)) <C E . Consequently, the significant contri- 
bution accumulates around the resonant energy level 



-En 



Eq — ^rf- 



(15) 



In terms of the generalised eigenstates, there will always 
be a continuous range of occupied states around the res- 
onant energy E Tes . 



III. VISIBILITY OF THE INTERFERENCE 
PATTERN 

The form of the free- falling atomic cloud \^u) was in 
Eq. (|12p expressed as a convolution of the outcoupling 
rate function f2(t) and a spatial term |$(i)). Next we 
will consider the emerging interference patterns due to 
multiple rf-fields driving the coupling simultaneously. 

The corresponding (classical) magnetic field compo- 
nents BJ f (t) interfere with each other, such that the total 
field is B r f(£) = J2i Brf(£). Correspondingly, the outcou- 
pling rate function f2(t) = ^ £li(t) and, because of the 
linearity of Eq. (|12[) , the outcoupled component is 



\*u(t)) = [o * (e|*»](t) - * (e|$»](t) 

i 

= £l*U*)>- (16) 

i 

The interference pattern appears similarly in the (quan- 
tum) matter fields as a sum of atomic streams, each of 
which corresponds to an atomic beam outcoupled by a 
single rf- field component. 



The point of view expressed in Ref. [33[ was that the 
combination of the (classical) magnetic fields, which op- 
erate at frequencies lo\ and W2 with equal constant am- 
plitudes, corresponds to a single field, whose carrier fre- 
quency is the average oj+ = (oji + 1^2) /2 and the pulse 
envelope is modulated at frequency oj_ = (loi —u)2)/2 (cf. 
Fig. [1]). Moreover, the relative phase difference between 
the circular components driving the outcoupling, 

A9 = arg(e Q ■ e^e^ 1 ) - arg(e Q ■ e^e 1 * 2 ), (17) 

shifts the envelope of the interference pattern and, con- 
sequently, the intensity profile of the falling stream of 
atoms. Generally, the interference pattern depends on 
(i) the relative amplitudes, (ii) the relative phase differ- 
ence, and (iii) frequency separation. 

Based on this description, one might expect that when- 
ever the carrier frequency u + falls into the region where 
the overlap integral \{il}E -hu + \4>)\ is finite, there would 
be a finite stream of atoms falling from the trap, and the 
intensity of the stream would be modulated at frequency 
u;_ , such that the maxima of the rf- field coincide with the 
maxima of the atomic intensity (see Fig. 3 in Ref. |33j). 



B. Interference of quantum fields 



The wave-packet result derived in Sec. Ill CI explains 
why the above-mentioned simplistic analogy from the 
classical interference is not exactly true. According to 
Eq. (TT5|) . the visibility of the interference pattern is af- 
fected by two contributions: (i) interference pattern of 
the magnetic fields and (ii) convolution by the tempo- 
ral free-fall evolution of the initial trapped state profile. 
Looking at the stream at a particular position x as a 
function of time, the interference pattern of the mag- 
netic fields, possibly with perfect visibility, is smoothed 
by the temporal width of the instantaneous outcoupled 
state |$(x, i)) falling past this point. 

For a Gaussian initial state \4>o), the analytical solution 



|$(x, t)) cx exp 



(x-x - hgt 2 ) 2 



2a{tf 



(18) 



where xq = g/oj 2 and a(t) = \/ (Tq + t 2 /a 2 , with co = 
yjh/muj, allows us to estimate the temporal width. 
Namely, at time t the wave packet has a spatial width of 
a(t) centralised around position xq + ^gt 2 , and the cen- 
ter of mass falls with velocity v(t) = gt, so the passing 
time is approximately a(t)/v(t) > l/crg; this value corre- 
sponds to a balance between dispersion and gravitational 
acceleration. Therefore, even with an infinitely long cou- 
pling time, the interference pattern is still smoothed by 
a distribution with a finite width. 

With a single rf-field, the amplitude of the falling 
atomic flux depends on the applied rf-frequency w r f . This 



TABLE I: The physical parameters used in the examples. 
Here, ao = 5.5 x 1CP 11 m is the Bohr radius. 



Quantity 


Symbol 


Value 


Trap frequency (x and z direction) 


0J x ,z/2lT 


160 Hz 


Trap frequency (y direction) 




6.7 Hz 


Rabi frequency 


|n|/27T 


50 Hz 


Bias frequency 


^bias/27T 


900 kHz 


Number of atoms 


N 


10 5 


Scattering length 


a 


103 a 



can be seen from the resonance energy condition (|15[) as 
compared to the overlap integral ("0b Therefore, two 
different rf-frequencies generally produce atomic streams 
with different amplitudes, even if the rf- field amplitudes 
are the same. According to Eq. P^)l . the relative phase 
difference of two magnetic fields (|17l) maps directly to the 
relative phase difference of the resulting matter waves. 
This was explicitly demonstrated in the experiment of 
Ref. [H. 



IV. APPLICATION AND COMPARISON TO 
NUMERICAL RESULTS 

In the following, we will concentrate on 87 Rb atoms 
and in particular its hyperfine ground state F = 1. In 
this case the Land factor is gF = —1/2 and, therefore, 
the trapped low-field-seeking Zeeman sublevel is |T) = 
\F = 1, Mp = —1) and the untrapped one is \U) = \F = 
l,Mp = 0) (see Fig. [T]). The physical parameters are 
adopted from Ref. [33j and are summarised in Table |TJ 

In this section we will compare the wave-packet solu- 
tion (|T12j) to numerical simulations including all the Zee- 
man sublevels. Especially, we show the impact of the 
atomic contact interactions by solving the corresponding 
Gross-Pitaevskii equation [4|, |5j . Both the model and the 
simulations are one dimensional. The contact interac- 
tions appear as an additional non-linear mean-field term 
Uint(x,t) — (x , t)\ 2 . The scaled interaction coeffi- 

cient is giu = (y/^i^2'rn/2'Kh)g3£) [iof . where the three- 
dimensional interaction term is <?3D = 4,nH 2 aN/m, with 
scattering length a and number of particles N 0, Q ■ 

When neglecting the atomic contact interactions, the 
ground state of the harmonic trapping potential is a 
Gaussian |(/> ). The overlap integral between the Gaus- 
sian state and the generalised energy eigenstates \i/je) 
can be calculated analytically [Til |23|. In the limit of 
a steep gravity slope, g 3> 0, the generalised energy 
eigenfunction approaches Dirac's delta distribution as 
\i/je(x)) ~ 5(x + E/mg)/y/mg. Since the width of the 
trapped state a clearly exceeds the characteristic length 
scale of the Airy distribution I, the overlap integral is well 
approximated by (V'_e|</ , o) — [ 7r ( m .9 cr ) 2 ]~ 1 ^ 4 exp[— (E + 
mgxo) 2 /2(mga) 2 }, as is obvious in Fig. [3J 

In Fig. [3] the time evolution of the instantaneous out- 
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FIG. 2: (Color online) The overlap integral {ipE -nu xt \<t>o), m 
arbitrary units, as a function of rf-frequency w r f. As men- 
tioned in the text, the form is well approximated by a Gaus- 
sian shape centralised at (E — mgxo)/h ~ 910.3 kHz with 
width mga/h ~ 1.8 kHz. 




1 1.5 2 

TIME (ms) 



FIG. 3: (Color online) The density profile of the instanta- 
neous outcoupled state |$(x,t)| 2 (arbitrary units). Initially 
the state is the Gaussian ground state of an harmonic po- 
tential. In the linear potential slope the functional form is 
maintained, while the center of mass accelerates according to 
classical mechanics, xo(t) 
perses, <r 2 (t) = <r 2 +£ 2 /cr 2 . 



x o + hgt 2 , and the width dis- 



coupled state (fTTj) for the non-interacting case is shown. 
The state corresponds to falling atoms outcoupled by an 
infinitesimally short rf-pulse. The total wave packet due 
to an rf-pulse with a finite duration, is achieved by inte- 
grating this state over time, in accordance with Eq. (|12l) . 
In the examples, we use a 5-ms-long box-shaped pulse 
form for the rf-field. The field amplitude is such that the 
maximum Rabi frequency is |OJ /2tt = 50 Hz. However, in 
the spirit of Ref. [33[ , we assume linear polarisation and 
the coupling is therefore suppressed by a factor l/v2- 
The density profile of the resulting stream of outcoupled 
atoms vs. time is shown in Fig. 2J The number of out- 
coupled atoms, as well as the density profile depends on 
the applied rf-frequency. This dependency is shown in 
Fig. [5] The increase in the density follows the amplitude 
of the overlap integral shown in Fig. [5J 

The density profile of a pulsating outcoupled stream, 
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FIG. 4: (Color online) Output-coupled atomic density for 
a 5-ms-long box-shaped pulse with rf- frequency cj r f/27r 
:! I u kHz. The density is in units of 10 3 1/m and is flattened 
from above to increase the clarity. 
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FIG. 6: (Color online) As Fig. [4] but with two simultaneous 
equally strong pulses with rf-frequencies u r f il /27r = 9f0 kHz 
and a>rf,i/27r =911 kHz, and with a relative phase difference 
of n. 
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FIG. 5: (Color online) Output-coupled atomic density profiles 
8 ms after the beginning of a 5-ms-long box-shaped pulse. The 
plot shows 8 different rf-frequencies; from lowest to highest 
density: w rf /27T = 907, 907.5, . . . , 910.5 kHz (cf. Fig. 0. 



FIG. 7: (Color online) As Fig. but with two simultaneous 
rf-fields with aj rf ,i/27r = 909 kHz and aj rf , 2 /27r = 908 kHz, 
and a relative phase difference of n. Analytical model (thick 
line) agrees well with numerical simulation (thin line). 



which is produced by two simultaneous resonant rf-pulses 
separated in frequency by Acj r f/27r = 1 kHz and in phase 
by A8 = it, is shown in Fig. |5] as a function of time. In 
Fig.[7]we compare the analytically calculated density pro- 
file to the numerically computed one. In the numerical 
computation the time-dependent Schrodinger equation 
was evolved taking into account three states, i.e., one har- 
monically trapped, one harmonically anti-trapped, and 
one affected by a linear potential with the slope corre- 
sponding to the gravity. The numerical computations 
were done for both interacting and non-interacting cases. 
Overall, we find good agreement between the analytical 
and the numerical results. 

Due to the atomic contact interactions the trapped 
ground state is broadened from a Gaussian into a 
Thomas-Fermi distribution. Accordingly, the range of 



rf-frequecies capable of producing outcoupling changes. 
In Fig. [5] we show how also in the interacting case the 
visibility of the interference pattern due to two equally 
strong rf-fields is not perfect. Especially, our example 
shows the interesting case, where one of the rf-frequencies 
outcouples hardly any stream while the other one and the 
average do. As in the non-interacting case, applying the 
equally strong fields simultaneously produces interference 
with low visibility. 

Finally, in Fig. [9j we compare the rf-pulses and in- 
duced outcoupling streams. The figure shows clearly how 
the visibility in the outcoupled atomic stream diminishes 
with increasing frequency separation in the causative out- 
coupling rf-fields, even if the rf-field itself is with perfect 
visibility. 
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FIG. 8: (Color online) As Fig. [S] but including the atomic 
contact interactions. Looking at the outcoupling near the 
edge of the distribution, there is a situation, where one rf- 
frequency (w r f,i/27r = 903 kHz, highest thin straight line) 
produces a strong stream of atoms, another rf-frequency 
(w r f : 2/27r = 901 kHz, lowest thin straight line) with the 
same field amplitude almost nothing, and the average rf- 
frequency (oj r f,i +w r f l 2)/2 (middle thin straight line) again a 
clear stream. Applying both equally strong fields simultane- 
ously, with a relative phase difference of ir, shows interference 
with a limited visibility (thick oscillating line) in accordance 
with the non-interacting examples. 
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FIG. 9: (Color online) The perfect interference patterns of 
the coupling magnetic fields (left panels) map to smoothed 
interference in the corresponding atomic density (right panels) 
because of the finite spatial extend of the trapped state. The 
used frequencies are w r f,i/27r = 911 kHz and (i) u I t,2/2ir = 
906 kHz (first row) (ii) u r f,2/27r = 908 kHz (second row), 
and (iii) cj r f,2/27r = 910 kHz (third row). The relative phase 
difference between the 5-ms-long pulses is n and the atomic 
density is plotted at time t = 8 ms. 



V. CONCLUSIONS AND DISCUSSION 



We have derived a linear wave-packet solution for out- 
put coupling scenario. The model establishes a bridge be- 
tween two different ways of looking at the interference of 
overlapping atom lasers, and shows that the effect can be 
understood equally as interference of spatially extended 
atomic clouds as well as interference of classical magnetic 
fields causing the output coupling. 

The model is built in terms of generalised energy eigen- 
states of a linear potential caused by gravity, and it shows 
how the total wave packet can be interpreted as being 
constructed by a continuous loading of a continuous spec- 
trum of these states, which individually do not corre- 
spond to a physical solution. In general, our model does 
not suffer from unphysical infinite quantities [U |3l| . 

Through the analysis of the solution, it was shown that 
the visibility of the observed interference pattern is lim- 
ited by the spatial extend of the trapped cloud, which 
serves as a source for the atomic beams. Furthermore, 
the visibility is shown to be affected by the rf-frequencies 
in the sense of selecting a resonant energy and, moreover, 
amplitude for the atomic stream. 

The simple linear model was then compared to numer- 
ical simulations including the atomic interactions as well 
as all the Zeeman sublevels, and the qualitative match 
was shown to be excellent using experimentally realistic 
parameters. The model is one dimensional and assumes 
weak coupling. The applicability is therefore restricted 
to cases, where the transversal extent of the source con- 
densate is wide 11, 23, 24]. Within these restrictions, the 
presented linear model also generalises straightforwardly 
for multiple dimensions and to an outcoupling scenario 
based on a Raman transition including an initial momen- 
tum kick. 
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